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Abstract 

We characterize non-degenerate Lagrangians of the form 

f(u x ,u y ,u t )dx dydt 

such that the corresponding Euler-Lagrange equations (f Ux ) x + (f Uy )y + (fu t )t = are integrable by 
the method of hydrodynamic reductions. The integrability conditions constitute an over-determined 
system of fourth order PDEs for the Lagrangian density /, which is in involution and possess 
interesting differential-geometric properties. The moduli space of integrable Lagrangians, factorized 
by the action of a natural equivalence group, is three-dimensional. Familiar examples include 
the dispersionless Kadomtsev-Petviashvili (dKP) and the Boyer-Finley Lagrangians, / = u x /3 + 
u 2 — u x u t and / = u 2 + u 2 — 2e Ut , respectively. A complete description of integrable cubic and 
quartic Lagrangians is obtained. Up to the equivalence transformations, the list of integrable cubic 
Lagrangians reduces to three examples, 

/ = u x u y u t , f = u x u y + UyU t and f — u x / 3 + u 2 - u x u t ( dKP ) . 

There exists a unique integrable quartic Lagrangian, 

f = u x + 2u 2 x u t - u x u y - u 2 . 

We conjecture that these examples exhaust the list of integrable polynomial Lagrangians which are 
essentially three-dimensional (it was verified that there exist no polynomial integrable Lagrangians 
of degree five). 

We prove that the Euler-Lagrange equations are integrable by the method of hydrodynamic 
reductions if and only if they possess a scalar pseudopotential playing the role of a dispersionless 
'Lax pair'. 
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1 Introduction 



The method of hydrodynamic reductions, originally proposed for the dKP equation in |12 | I13 | HI], 
applies to a broad class of multi-dimensional dispersionless PDEs F(u, u s , . . .) = where u is a (vector- 
) function of d independent variables s = (x,y,t, . . .). The main idea is to look for solutions in the 
form u = u(R 1 , . . . , R n ) where the Riemann invariants R 1 (x, y, t, . . .), . . . , R n (x, y,t, . . .) are arbitrary 
solutions of (d — 1) commuting diagonal systems 

Ry = fJ, l (R) R l x , R\ = \ % {R) R l x , 

etc. We recall, see |31] . that the commutativity conditions imply the following restrictions on the 
characteristic speeds: 



Thus, the original multi-dimensional equation is decoupled into a collection of commuting n-component 
(1 + l)-dimensional systems in Riemann invariants which can then be solved by the generalized hodo- 
graph method |31| . Solutions arising within this approach, known as nonlinear interactions of n planar 
simple waves, were extensively investigated in gas dynamics (simple waves, double waves, etc, |29] ) 
and the theory of dispersionless (2 + l)-dimensional hierarchies of the dKP type El EJ El I2T1 

EH ESEOl 0123 CHI EH- 

It was suggested in |U] to call a multi-dimensional equation integrable if it possesses 'sufficiently 
many' n-component hydrodynamic reductions (parametrized by (d — 2)n arbitrary functions of a single 
variable). It turned out that this requirement is very strong and provides the effective classification 
criterion. Partial classification results for (2 + l)-dimensional integrable systems of hydrodynamic 
type, 

u t + A(u)u x + B(u)u y = 0, 

(here u is m-component vector, A and B are mxm matrices) were obtained in El - It was observed 
that the class of PDEs amenable to the method of hydrodynamic reductions extends beyond the class 
of hydrodynamic type systems. For instance, the method applies to scalar second order PDEs 

Fyllxxi ^xyt ^xti ^yyi l^yti ^Ht) — 0j 

see e.g. j2H]; HHJ for the classification of integrable equations of the form uu = f(u xx ,u x t,u xy ). 
Hydrodynamic reductions of multi-dimensional dispersionless integrable systems in dimensions greater 
then three were discussed in (HI El • 

In this paper we apply the method of hydrodynamic reductions to three-dimensional equations of 
the form 

(fuj x + (fu y )y + (fu t )t = 0, (2) 

which are the Euler-Lagrange equations corresponding to first order Lagrangian densities f(u x ,u y , ut). 
We begin with two illustrative examples. 

Example 1. Consider the linear wave equation Uu = u xx + u yy corresponding to the quadratic 
Lagrangian density / = u x + u y — . Setting a = u x , b = u y , c = ut we can rewrite it in the first 
order form 

a y = b x , at = c x , b t = c y , a x + b y - c t = 0. (3) 

Let us seek solutions in the form a = a(R . . . , R n ), b = b(R 1 , . . . , R n ), c = c{R\ . . . , R n ) where 
the Riemann invariants R 1 (x,y,t), R n (x,y,i) are arbitrary solutions of a pair of commuting 
hydrodynamic type systems 

R y = fi l (R) Rl, R x = \ % {R) R\\ 
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(for the sake of symmetry we have chosen t as a distinguished variable) . Substituting this ansatz into 
© one obtains the equations dib = fi l diC, did = X l diC along with the dispersion relation 

(A*) 2 + Ox*) 2 = 1. 

Setting A* = cos ip l , // = sin ip % and taking into account the commutativity conditions (^Q) one obtains 
djip 1 = for any j ^ i. Thus, ip l = ip l (R l ) so that both systems of hydrodynamic type assume 
decoupled forms 

R l y = sin ^ (if) R\, R x = cos ^ (IV) R\, 
with the general solution R l (x,y,t) given by the implicit relation 

g^R 1 ) = t + ysmtp i (R i ) + xcos<p i {R i ); (4) 

here g l (R l ) are yet another n arbitrary functions. Finally, the consistency conditions for the equations 
dib = fi l diC, did = X l diC imply didjC = so that one can set c = R 1 + ... + R n (notice that one has a 
reparametrization freedom R % — * f l {R 1 )). Thus, we arrive at solutions of the wave equation given by 
the formulae 

ut = c = R 1 + ... + R n , 
u y = b = f smip 1 (R l )dR 1 + ... + / smip n (R n )dR n , 

u x = a = f cos tp\R x )dR x +...+/ cos <p n (R n )dR n , 

where R l (x,y,t) is defined by the implicit relation These solutions depend on 2n arbitrary 

functions of a single argument and can be viewed as linear superpositions of n elementary solutions 
of the form 

ut = R, u y = j sia.<f(R)dR, u x = f cos (p(R)dR 

where R(x,y,t) is defined by the implicit relation g(R) = t + y sin. <p(R) + x cos (p(R). This relation 
implies that the level sets R = const are null-planes so that the gradient of u is constant along a 
one-parameter family of null-planes. Solutions of this type are known as planar simple waves. 

Example 2. Let us apply the same approach to the nonlinear wave equation e Ut uu = u xx + u yy 
(also known as the Boyer-Finley equation) which corresponds to the Lagrangian density / = u x + Uy — 
2e Ut . Setting a = u x , b = u y , c = ut we can rewrite it as follows: 

a y = b x , a t = c x , b t = c y , a x + b y - e c c t = 0. (5) 

We again seek solutions in the form a = a(R , . . . , R n ), b = b(R x , . . . , R n ), c = c(R l , . . . , R n ) where 
the Riemann invariants solve hydrodynamic type systems 

R y = n\R) R\, R l x = \ % {R) R\- 

The substitution into © implies 

dib = //<9jC, did = X l diC (6) 

along with the dispersion relation 

(V) 2 + Gt/) 2 = e c . 

Setting X 1 = e c l 2 cos <*p l , fj? = e C//2 sin (p l and taking into account the commutativity conditions Q one 
obtains the expression for djip 1 : 

d ^ = A cot t_^ d ^ i^j. (7) 
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The consistency conditions for the equations (JBJ) imply 



didjC = —r^ — r^diC + —r^ — r^djC, i ^ j. 
\o - A 1 A 1 - \i J ^ J 

Substituting here X 1 = e C//2 cos p l and taking into account ((7| we obtain the following system which 
governs hydrodynamic reductions of the nonlinear wave equation: 



1 ip l -ip> djcdic 
dup = — cot - — djc, djdic = - — : — -: (8) 

2 2 2 sin 2 y y 



notice that this system is a nonlinear analog of the system dj(f l = 0, djdiC = derived in Example 
1 for the linear wave equation. One can show that the system (|SJ) is in involution (that is, the 
compatibility conditions dk(dj(p l ) = dj(dkf l ) and dk(djdic) = dj(dkdic) are satisfied identically) and 
its general solution depends, modulo reparametrizations R 1 —* f l (R l ), on n arbitrary functions of a 
single argument. Once a solution to the system (jHJ) is found, one can reconstruct b{R) and a(R) from 
the equations (JHJ) which are consistent by construction. After that one has to solve hydrodynamic 
type systems (where the characteristic speeds ^ 1 {R) and A 1 (it!) are known) which can be done by 
the generalized hodograph method [31]. This gives some more n arbitrary functions. Thus, solutions 
arising from n-component reductions of the nonlinear wave equation depend on 2n arbitrary functions 
of a single argument. For one-component reductions, equations (JHJ become vacuous and, setting c = R, 
one has 

u t = c = R, u y = b = J e R l 2 sin<p(R)dR, u x = a = f e R l 2 cos Lp(R)dR 

where the function R(x, y, t) is defined by the implicit relation g(R) = t+ye R / 2 simp(R)+xe R ^ 2 cos <p(R)- 
This relation implies that the level sets R = const are characteristic planes. Solutions of this type 
are known as planar simple waves. Thus, solutions governed by the system (jHJ) can be interpreted as 
nonlinear interactions of planar simple waves. We refer to H3 1221 119j for further discussion and 
explicit examples. 

The main object of our study are nonlinear Lagrangian PDEs (of the type discussed in Example 
2) which are integrable in the above sense, that is, possess infinitely many hydrodynamic reductions. 
This requirement imposes strong restrictions on the Lagrangian density /. Our first main result is the 
system of partial differential equations for the density f(a,b,c) providing the necessary and sufficient 
conditions for the integr ability (we use the notation a = u x , b = u y , c = u t ). These conditions can 
be represented in a remarkable compact form: 

Theorem 1 For non- degenerate Lagrangians, the Euler- Lagrange equations HJ) are integrable by the 
method of hydrodynamic reductions if and only if the density f satisfies the identity 

dJ-f S 

d A f = d 3 f— + -det(dM). (9) 

Here d 3 f and <i 4 / are the symmetric differentials of / which appear in the standard Taylor expansion 
f(s + ds) - f(s) = df + d 2 f/2\ + d 3 //3! + ... for a function f(s) of three variables s = (a, b, c). The 
Hessian H and the 4x4 matrix M are defined as follows: 



faa fab fa* 

U d< I j f ab f bb f bc I . .1/ 

fac fbc fa 



( 


fa 


fb 


fc \ 


fa 


faa 


fab 


fac 


fb 


fab 


fbb 


fbc 


\ fc 


fac 


fbc 


fee j 



(10) 



4 



The differential dM = M a da + M^db + M c dc is a matrix-valued linear form 



/ n 
\ 


f 

J aa 


f , 

J ab 


f 

J ac 


\ 
] 




/ n 
I 


f , 
J ab 


Jbb 


f, 

Jbc 


\ 
] 




/ n 


f 

J ac 


f, 

Jbc 


f 

J cc 


\ 
\ 


J aa 


f 

J aaa 


f u 
J aab 


f 

J aac 




da + 


J ab 


f u 
J aab 


J abb 


J abc 




db + 


J ac 


f 

J aac 


J abc 


f 

J acc 




fab 


faab 


fabb 


fabc 




fbb 


fabb 


fbbb 


fbbc 




fbc 


fabc 


fbbc 


fbec 




\ fac 


faac 


fabc 


face 


J 




\ fbc 


fabc 


fbbc 


fbec 


J 




\ fee 


face 


fbec 


fece 


) 



dc. 



A Lagrangian is said to be non-degenerate iff H ^ 0. 

Both sides of the equation are homogeneous symmetric quartic forms in commuting formal 
variables da, db, dc. Equating similar terms we obtain expressions for all fourth order partial derivatives 
of / in terms of its second and third order derivatives (15 equations altogether). The resulting over- 
determined system for / is in involution, and its solution space is 20-dimensional (indeed, the values 
of the derivatives of / up to order 3 at a given point (ao,6(b c o) amount to 20 arbitrary constants). 
Factorised by the action of a natural equivalence group of dimension 17 (generated by arbitrary affine 
transformations of the independent variables a,b,c plus transformations / — > [if + aa + fib+'jc+S, see 
the end of Sect. 2 for the discussion of the origin of these symmetries) this provides a three-dimensional 
moduli space of integrable Lagrangians. Details of the derivation of the integrability conditions @ are 
given in Sect. 2. Notice that these conditions arise in a somewhat different form, namely, as explicit 
formulae for the fourth order derivatives of / (which are very complicated). It is a truly remarkable 
fact that they compactify into a single expression (jHJ). Differential-geometric aspects and an invariant 
formulation of the integrability conditions (PJ) are discussed in Sect. 5. 

Example 3. For the dKP Lagrangian (/ = a 3 /3+b 2 -ac) one has d A f = 0, H = -2, det(dM) = 0. 
Similarly, for the Boyer-Finley Lagrangian (/ = a 2 + b 2 — 2e c ) one has d 4 f = — 2e c (dc) 4: , d 3 f = 
—2e c (dc) 3 , H = — 2e c , det(dM) = 0. In both cases the identity © is obviously satisfied. 

Remark. Notice that in two dimensions any Euler-Lagrange equation of the form (fu x )x+(fu y )y 
is automatically integrable. Indeed, in the new variables a = u x , b = u y it takes the form of a 
two-component quasilinear system a y = b x , (f a )x + (fb)y = which linearises under the hodograph 
transformation interchanging dependent and independent variables. This trick, however, does not 
work in more than two dimensions. 

Sect. 3 is devoted to polynomial Lagrangians. Our first result is a simple Lemma stating that, 
for non-degenerate homogeneous solutions of the equation Q, the degree of homogeneity can take 
only one of the three values 0, 2 or 3 (recall that Lagrangians of homogeneity one are automatically 
degenerate). This observation is particularly useful for the classification of polynomial Lagrangians 
implying that the 'leading' homogeneous part of a polynomial solution / must be either of degree 3 
or degenerate (polynomials of degree two give rise to linear Euler-Lagrange equations). In Sect. 3.1 
we obtain a complete list of integrable cubic Lagrangians with the densities 

/ = C(u x , u y , u t ) + au 2 + f3u 2 + 'juf + ijlu x u v + vu x u t + r/u y u t ; 

here C is a homogeneous cubic form. The substitution into the integrability conditions © implies 
that C must be totally reducible, leading to three essentially different possibilities: 

(1) C consists of three lines in a general position. Further analysis allows one to eliminate quadratic 
terms, leading to a unique Lagrangian density 

/ = u x u y u t 

with the corresponding Euler-Lagrange equation UfU xy + u y u x t + u x u y t = 0. Although this example 
looks deceptively simple, the corresponding dispersionless Lax pair is quite non-trivial, see formula 
(|12|) . To the best of our knowledge this Lagrangian has not been recorded before. 

(2) C contains a double line. Up to equivalence transformations this case reduces to the density 

/ = u 2 u y + UyUt 
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which generates a particular flow of the so-called r-Dym hierarchy. 
(3) C is a triple line. This case reduces to the dKP Lagrangian density 

/ = ul/3 + u 2 - u x u t . 



Quartic Lagrangians are classified in Sect. 3.2. Here the leading part of the density / must 
necessarily be degenerate, leading to a unique integrable example 

/ = ?4 + 2u 2 x u t - u x u y - uf 

which corresponds to a particular flow of the so-called r-dKP hierarchy. 

We have verified that there exist no polynomial solutions to the system © of degree five. It is 
tempting to conjecture that the four examples listed above exhaust the list of polynomial integrable 
Lagrangians. 

Homogeneous Lagrangians can be obtained by setting f(a,b,c) = a k g(^,rj), £ = b/a, rj = c/a 
(recall that the degree of homogeneity k can take the values 0, 2, 3 only). Substituting this ansatz into 
the integr ability conditions @ and analysing the resulting equations for g(^,rj) one can show that the 
case k = 3 leads to the cubic Lagrangian / = u x u y ut and, thus, gives no new examples. The case 
k = 2 reduces to functions g which are quadratic in £, rj and, hence, generate quadratic Lagrangians 
with linear Euler-Lagrange equations. The last case k = turned out to be quite nontrivial. The 
substitution of f(a,b,c) = g(£,rj) into (jHJ) results in a system of five equations expressing all fourth 
order partial derivatives of g in terms of lower order derivatives. In symbolic form, this system can be 
represented as follows: 

d 4 g = d 3 g^ + Aet(dm) + 3^-det(dn). (11) 
h h h 

Here, as in @, d s g are symmetric differentials of g (notice that g is now a function of two variables), 
the matrices m and n are defined as 



g$ g v 
I" | 9i 9tt 9^ 

9r\ 9€j\ 9r\r] 



9£,ri 9t)r] 



and 

h = -det(m) = g 2 g^ - 2g^g v g^ + g^g^- 

The non-degeneracy of the Lagrangian density f(a, b, c) = g(£, rj) is equivalent to the condition h ^ 0. 
One can show that the over-determined system (jlljl is in involution and its solution space is 10- 
dimensional (indeed, the values of partial derivatives of g up to order 3 at a fixed point (£o>%) 
amount to 10 arbitrary constants). Although it is still difficult to integrate this system in general, 
some particular solutions can readily be constructed. For instance, any quadratic form g(£, rj) = 
a£ 2 + P^rj + 77/ 2 + jti£ + vrj is a solution of (|11|) , leading to integrable Lagrangians with the densities 

U 2 Ut Uy U 2 Uy U t 

f = a-f +13 — f- + 7-| + + v — • 

'^x 

Up to the equivalence transformations any expression of this type can be reduced to either of the 
non-equivalent (over reals) canonical forms, 

y} y + uf _ UyU t _ «| ut_ 

J 2 ' ^ 2' * 2 
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In Sect. 4 we study scalar pseudopotentials for integrable Euler-Lagrange equations (also known 
as S- functions, dispersionless Lax pairs, etc). Examples thereof include the dispersionless Lax pair 



St — ^ S x + u x S x + iiy, Sy — ^ S x + u x 



which generates the dKP equation u x t — u x u xx = u yy . Similarly, the Boyer-Finley equation u xy = {e Ut )t 
possesses the dispersionless Lax pair 



S t = u t - In S y , S x = u x 



Sy 



Further examples of integrable multi-dimensional equations possessing pseudopotentials of the above 
type can be found in |32| \2H \ IT7| . It was proved in [H] that for two-component (2 + l)-dimensional 
systems of hydrodynamic type the existence of dispersionless Lax pairs is necessary and sufficient for 
the integrability (that is, for the existence of sufficiently many hydrodynamic reductions). Dispersion- 
less Lax pairs constitute a basis for the twistor and dispersionless <9-approaches to multi-dimensional 
dispersionless hierarchies (2H1 El I2j • 

Our second main result is the following 

Theorem 2 The Euler-Lagrange equation (f Ux )x + (fu y )y + (fu t )t = is integrable by the method of 
hydrodynamic reductions if and only if it possesses a dispersionless Lax pair 

S t = F(S X , U x , Uy, u t ), Sy = G(S X , u x , Uy, u t ). 

In some cases it seems to be more convenient to work with parametric Lax pairs 

St = F(p, u x , u y , u t ), S y = G(p, u x , u y , u t ) , S x = H(p, u x , u y , ut), 

which take the above form if one expresses the parameter p in terms of S x from the third equation. For 
instance, the equation Utu xy + u y u x t + u x u y t = corresponding to the Lagrangian density / = u x u y ut 
possesses the parametric Lax pair 

— = C(p)» — = C(p) + r~\ ' — = C(p) + o r \ ; ( 12 ) 

u x Uy 2p(p) u t 2p(p) 

here (p') 2 = 4p 3 + A 2 and (' = —p (Weierstrass p and ( functions, see Sect. 3). 



Differential-geometric aspects of the integrability conditions © are investigated in Sect. 5. The 
main object associated with the Lagrangian density /(a, b, c) is the Hessian metric d 2 f. We show that, 
by virtue of @, this metric is necessarily conformally flat. 

Remark. Notice that any Euler-Lagrange equation (f Ux )x + (fu y )y + (fu t )t = is manifestly 
conservative and, moreover, possesses three extra conservation laws 

(U x fu x - f)x + (Uxfu y )y + (u x f Ut ) t = 0, 



(Uyfu x )x + (u y fu y ~ f)y + {Uyfu t )t = 0, 
{Utfujx + (Utfu y )y + {Utfu t ~ f)t = 0, 
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which are components of the energy-momentum tensor. In the dKP case, / = u x /3 + u 2 — u x ut, these 
four conservation laws take the form 

(Ul - Ut) X + (2Uy)y - (U X ) t = 0, 

(2u x /3 - u 2 y ) x + {2u x u y )y - (u 2 x ) t = 0, 

{u 2 u y - u y u t ) x + (u x iit + u 2 - u x /3) y - (u x u y ) t = 0, 

(u 2 x u t - u 2 ) x + (2u y u t )y - {v 2 + «3/3) t = 0, 

respectively. One can show that the dKP equation possesses no extra conservation laws of the form 

g(u x ,u y , u t )x + h(u x , u y , u t ) y + p(u x , u y , u t ) t = 0. 

The discussion of the corresponding hierarchies of higher nonlocal symmetries and conservation laws 
for integrable Euler-Lagrange equations is beyond the scope of this paper. Our primary goal is the 
characterization of integrable Lagrangians based on the method of hydrodynamic reductions. 

2 Derivation of the integrability conditions: proof of Theorem 1 

Introducing the variables a = u x , b = u y , c = ut one rewrites the Euler-Lagrange equation © in the 
first order form 

a y = b x , a t = c x , b t = c y , (f a ) x + (f b ) y + (f c ) t = 0. (13) 

The idea of the method of hydrodynamic reductions is to look for solutions of the system (|13[) in 
the form a = a(R , . . . , R n ), b = b(R \ . . . , R n ), c = c(R , . . . , R n ) where the Riemann invariants 
R 1 (x,y,t), . . . , R n (x,y,t) are arbitrary solutions of a pair of commuting hydrodynamic type flows 

R y = n\R)Rl R\ = \ i {R)R i x . (14) 

Substituting this ansatz into lfT5|) one obtains the equations 

d{b = fj. l d{a, die = X l dia (15) 

(here di = d/dpi) along with the dispersion relation 

L>(A\ //) = f aa + 2f abt S + 2f ac \ l + W) 2 + 2f bc n i \ i + f cc {\ 1 ) 2 = 0. (16) 

Hereafter, we assume the conic (|16[) to be irreducible. This condition is equivalent to the non- 
degeneracy of the Lagrangian or, equivalently, to the non- vanishing of the Hessian: H ^ 0, see (|10JI. 
The consistency conditions of the equations Q15|) imply 

<9'A* d\3 
didja = — ^ — r^dia + — -rdta. (17) 

Differentiating the dispersion relation ()16j) with respect to j ^ i, and keeping in mind (|15|) and 
((T)l one obtains the explicit expressions for djX" and djfi 1 in the form 

djX = (A* - \ J )B, ,<),„. d jf / = {fi l - n j )Bijd ja (18) 

where Bij are rational expressions in X", X J , //, fx 3 whose coefficients depend on partial derivatives of 
the density f(a,b,c) up to third order. Explicitly, one has 

B = % = ^ . 

13 D i3 2(f aa + + 0) + f^X* + X) + hbn^i + fbc(^X j + V 3 X) + f cc X V) ' 
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notice that, modulo the dispersion relation (|lfi)L the denominator Dij equals 4D y A ^ , i^ -4p-j • The 
numerator Nij is a polynomial expression of the form 

N i3 = f aaa + fa^ifjP + 2//) + f aac (\i + 2A*) + /«»mV + + /accA^A* + 2A')+ 
2/„& c (AV + AV + AV) + fbbbW) 2 ^ + fccc(\ l ) 2 \ J + 
/ 66c //(AV + 2AV) + //,,, A'lA'/r' + 2AV)- 
The equation (|T7|) takes the form 

(9j0ja = — {Bij + Bji)diadja. (19) 

The compatibility conditions dkdjX 1 = djd k X l , d k djfi l = djd k fi l and dkdjdia = djd k did are equivalent 
to the equations 

dfe-By = (BijBkj + BijBik — B k jB ik )d k a, (20) 

which must be satisfied identically by virtue of (|15|) . (jlfij) . The details of this calculation (which 

is computationally intense) can be summarized as follows. 

In order to obtain equations with 'smallest possible' coefficients at the fourth order derivatives of 
f(a,b,c) we rewrite (|2U)l as 

d k Nij = (jj-dkDij + B kj d k a + B ik d k a^j - DijB kj B ik d k a. (21) 

The fourth order derivatives of f(a,b,c) are present only in the l.h.s. term d k Nij. Further reduction 
of the complexity of the expression in the r.h.s. is achieved by representing 1/Dij in the form 

= Uij = [2(AV6c + fab)(Xfbc + fab) ~ fbb^^fcc + (A* + X j )fac + faa) 
+fbb(X J fbc + fabW + MX'fbc + fab)^ 3 + / fe >V]/((A 4 - X j ) 2 H) 

(which holds identically modulo the dispersion relation (|16|l). and the subsequent substitution B st = 
N s t/D s t = N s tU s t- The denominators of the r.h.s. terms in (|21jl cancel out as explained in the 
attached program file 2-Lagr3dim.frm producing a polynomial in A 1 , X J , X h , //, n J , fi k with coefficients 
depending on the derivatives of the density /(a, b, c). This was the most essential technical part of the 
calculation: the starting expression for the r.h.s. of ()21j) has more than 500.000 terms with different 
denominators; after properly organized cancellations we get a polynomial expression with less than 
2000 terms. Using (|16|) and assuming f^b ^ (this can always be achieved by a linear change of 
independent variables x,y,t), we simplify this polynomial by substituting the powers of (fi l ) s , (^Y, 
(l^ k ) s , s > 2, arriving at a polynomial of degree one in each of //, ^ 1 , [i k and degree four in A's. Equating 
similar coefficients of these polynomials in both sides of (|21|) we arrive at a set of 45 equations for 
the derivatives of the Lagrangian density f(a,b,c) (linear in the fourth derivatives). Solving it, we 
get closed form expressions for all fourth order derivatives of f(a, b, c) in terms of its second and third 
order derivatives. A straigtforward computation shows that the compatibility conditions are satisfied 
identically. Although the arising expressions are very long indeed, they can be rewritten in a compact 
form 0. This finishes the proof of Theorem 1. 

Further particulars and the programs in FORM |33j and Maple 1 [HI] for the computations described 
above are given in the attached files. The file 1-README explains the overall structure of this program 
package. 

Maple(TM) is a trademark of Waterloo Maple Inc. 
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Remark. The integrability conditions © are invariant under the obvious equivalence transfor- 
mations generated by 

(el) linear transformations of a, b, c corresponding to linear changes of the independent variables x, y, t; 

(e2) translations in a, b, c corresponding to the transformation u — * u + ax + (3y + jt; 

(e3) transformations of the form / — * fif + va + rjb + rc + p which do not effect the Euler-Lagrange 

equations. 

These transformations generate a 17-parameter symmetry group of the system 0. Notice that 
(el) and (e2) generate the group of affine transformations. These symmetries allow one to considerably 
simplify the classification results. For instance, given any quadratic polynomial Q{a, b, c) and a linear 
function l(a,b,c), both not necessarily homogeneous, the density 

f(a, b, c) = Q(a, b, c)/l(a, b, c) 

is a solution of Q- Up to the equivalence (el) - (e3) this expression can be transformed to either of 
the non-equivalent (over reals) canonical forms 

b 2 + c 2 + 1 be + 1 b 2 + c 

3 3 * 

a a a 

The classification below is carried out up to this natural equivalence. 



3 lntegrable polynomial Lagrangians 



We will start with a useful remark on homogeneous solutions of the equation (jJJJ • Recall that a function 
/ is said to be homogeneous of degree k if it satisfies the Euler identity af a + bfb + cf c = kf. This 
implies, in particular, that f a , ft,, f c are homogeneous of degree k—1, that is, 




(fc-1) 




(22) 



here the 3x3 matrix G is the Hessian matrix of /. 

Lemma Let f be a homogeneous solution of the equation with non-zero Hessian H = det G. 
Then the degree of homogeneity k can take only one of the three values 0, 2 or 3. 

Proof: 

The equation ® is an identity in da, db, dc. Let us replace da, db, dc by a, b, c, respectively. Under 
this identification one has: df -> kf, d 2 f (k - l)kf, d 3 f — ► (k — 2)(k — l)kf, d A f —>■(& — 3)(k - 
2)(k-l)kf. Moreover, 



dH 



H a da + H b db + H c dc H a a + H h b + H c c 



3{k 



H H 

since H is homogeneous of degree 3(k — 2). Finally, the matrix one-form dM reduces to 

/ (k-l)f a (k-l)f b (k-l)f c \ 

(k-l)f a (k-2)f aa (k-2)f ab (k-2)f ac 
(k-l)f b (k-2)f ab (k-2)f bb (k-2)f bc 
\(k- l)/ c (k - 2)f ac (k - 2)f bc (k - 2)f cc J 

its determinant equals (k — l) 2 (k — 2) 2 det M. Thus, for homogeneous /, the equation © implies 



(k - 3)(Jb - 2){k - l)kf = 3(k - 2) 2 {k - l)kf + 



3(k-l) 2 (k 
H 



-det M. 



(23) 
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We also point out the identity 

det M 
H 

Taking into account (|22|). this implies 

(fc-l)^ = -(/. A / c )| & \=-kf- 

Thus, both terms on the right hand side of (j23j) cancel, leaving the identity (k — 3)(k — 2)(k — l)k = 0. 
It remains to point out that Lagrangians of homogeneity k = 1 are automatically degenerate. Q.E.D. 

This result is particularly useful for the classification of polynomial solutions. Indeed, let us 
represent a polynomial / of degree k in the form 

f = Qk + Qk-i + ■■■ 

where Qj~ is a homogeneous polynomial of degree k in the variables a, b, c. Writing the equation Q in 
the homogeneous form 

Hd 4 f = d 3 fdH + 3det(dM), 

we readily see that the leading term Q/% must be a solution itself. This leads to two possibilities: either 
k = 3 (see Sect. 3.1) or, if k > 4, the Hessian of Qt must vanish identically. In the last case the 
classical result ^S] (see also [2f)], p. 234) implies the existence of a linear change of variables a,b,c 
after which Qk becomes a binary form, that is, a homogeneous function of two variables only (say, a 
and b). This partial 'separation of variables' considerably simplifies all calculations (see Sect. 3.2 for 
a complete analysis of the case k = 4). 



3.1 Classification of cubic Lagrangians 

The integrability conditions © provide a straightforward classification of integrable cubic Lagrangians 
with the densities 

f(a, b, c) = C(a, b, c) + era 2 + [3b 2 + 7c 2 + fiab + vac + rjbc 

where C is a homogeneous cubic form in a, b, c. Using equivalence transformations of the type (el) 
one can bring C to a canonical form, thus simplifying the analysis. Direct calculations using Q reveal 
that the cubic C must necessarily be totally reducible, that is, a product of three linear forms (this 
condition is independent of the quadratic part). Thus, we have four cases to consider, depending on 
the mutual position of the corresponding three lines. 

Case 1: Three lines in a general position. Without any loss of generality we can assume C = abc. 
The substitution of / into the integrability conditions Q readily implies a = (3 = 7 = 0. Since the 
remaining constants fJ>,v,rj can be eliminated by the equivalence transformations (e2), we have just 
one integrable Lagrangian density in this class, 

/ = U X UyU t , 

with the corresponding Euler-Lagrange equation 

U t U X y + UyU xt + U X Uy t = 0. (24) 
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Notice that this Lagrangian density is equivalent to 

/ = u 3 x + u 3 + u 3 - 3u x u y u t ; 

the corresponding (complex) linear change of variables x, y, t can be reconstructed from the factoriza- 
tion 

1 \/3 

a 3 + b 3 + c 3 - 3abc = (a + b + c)(a + eb + ec)(a + e& + ec), e = -- + i— ; 

in this case two of the three lines are complex conjugate. One can show that the equation (|2*i|) possesses 
the Lax pair 

S t = q{S x /u x ) u t , S y = r(S x /u x ) u y , (25) 
where the functions q(s) and r(s) (s = S x /u x ) satisfy a pair of ODEs 

Q—r , r - q 



r 



This system can be integrated in elliptic functions as follows. Introducing f(s) = q(s) — s and 
g(s) = r(s) — s one first rewrites the system in a translation-invariant form 

f' = f/g-2, g' = g/f-2 

or, equivalently, 

(fg)' = -(f + g), (/ + s)' = -^4^-- -6. 

fg 

Setting fg = u,f + g= —u' one obtains a second order ODE for u, namely, uu" + (u 1 ) 2 — 6u = 0, 
whose integral is given in parametric form u = p(p), s = ((p) (that is, u(s) is obtained by excluding p 
from these two equations). Here p(p) is the Weierstrass p-function, (p'(p)) 2 = 4p 3 (p) + A 2 (notice that 
92 = 0, 53 = —A 2 ), and ((p) is the corresponding zeta-function: C{p) = ~p(p)- Thus, in parametric 
form, fg = p(p), f + g = p'(p)/p(p)- The Lax pair (|23|) can be rewritten in parametric form as 
follows: adding and multiplying equations l(23|) we have 

^ + ^=q + r = 2s + f + g = 2ap) + ^f\, 

U t Uy P[P) 

^ = qr = s 2 + s{f + g)+fg = C 2 ^ + C (p)^M + 

U t Uy P(P) 



Solving these equations for ^ and ^ and keeping in mind that ^ = s = C(p)> we arrive at parametric 
equations 



u x u y 2p(p) u t 2p(p) 

Notice that these equations imply the algebraic identity 

U x UyJ \U X Uf 

Case 2: Three lines through a common point. Without any loss of generality we can assume 
C = ab(a + b). The substitution of the corresponding / into the integrability conditions Q implies 
7 = ^ = 77 = so that there is no c-dependence. Thus, this case gives no non-trivial examples. 
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Case 3: One double line. Here C = a 2 b, and the substitution into the integrability conditions 
implies 7 = v = 0. Translations in b and a eliminate a and /i. Furthemore, the linear transformation 
Pb + rjc — > c reduces / to the canonical form / = bc + a 2 b. The corresponding Lagrangian density 

/ = u 2 x u y + u y u t 

generates the equation 

u yt + u y u xx + 2u x u xy = (27) 
which, up to a rescaling, is a particular form of the so-called dispersionless r-Dym equation pQ, |2(J| . 

3 - r / 1 1 



^irf n In Uyltxx -. U x 11 x y 

y 2 - r \2 - r 9 1 — r 

Indeed, for r = 4/3 this equation possesses the Lagrangian 

/15 
(u y u t — —u^Uy) dxdy. 

The equation (JT7j) possesses the Lax pair 

St = -2u x S x + ^5^/2, S y = 2u y S x 1 / 2 . 

Case 4: One triple line. Here C = o 3 /3, and the substitution into the integrability conditions 
implies a single constraint 4/?7 — rj 2 = 0. Eliminating a by a translation of a, one arrives at the 
expression / = a 3 /3 + {yffib + ^/7c) 2 + a(fj>b + i^c). The linear change of variables \ffib + ^/jc — ► 
6, [lb + vc —> —c reduces / to the canonical form / = a 3 /3 + 6 2 — ac which corresponds to the dKP 
density 

/ = u 3 x /3 + u 2 - u x u t . 

Geometrically, the cases 3 and 4 can be viewed as degenerations of the case 1. It would be 
interesting to perform these degenerations explicitly on the level of the corresponding PDEs and Lax 
pairs. 

3.2 Classification of fourth order Lagrangians 

According to the Lemma, the leading quartic part of the Lagrangian density / must necessarily be 
degenerate and, thus, can be written as a form in two variables (say, a and b) after an appropriate 
linear transformation. Since any homogeneous binary quartic can be reduced to one of the five non- 
equivalent forms a 4 + fia 2 b 2 + 6 4 , a 2 b(a + b), a 2 b 2 , a 3 6 or a 4 (depending on the mutual location of 
its four roots), we have five cases to consider. A direct substitution into @ implies that the first 
four cases lead to inconsistency (whatever cubic and quadratic terms are), while the last case leads, 
up to the equivalence transformations (el)-(e3), to a unique solution / = a 4 + 2a 2 c — ab — c 2 . The 
corresponding Lagrangian density 

/ = u\ + 2u 2 u t — u x u y — u 2 

generates the equation 

u xy = -u t t + 2u t u xx + 4u x u xt + §u 2 u xx . 
Up to a rescaling, this is a particular case of the so-called r-th dispersionless modified KP equation 

3-r (3-r)(l-r) (3 - r)r (3-r)(l-r) 2 

U>xy — 77t \o^Ht n Ut'U'xx ^x^xt ^x^xxt 

(2 — r) z 2 — r 2 — r 2 

corresponding to the parameter value r = 2/3 (for which the equation becomes manifestly Lagrangian). 

A similar analysis reveals that the integrability conditions @ possess no polynomial solutions of 
degree 5. We conjecture that the four examples listed above exhaust the list of polynomial integrable 
Lagrangians which are essentially three-dimensional. 
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4 Pseudopotentials: proof of Theorem 2 



In this section we prove that any integrable non-degenerate Euler-Lagrange equation (J2J) possesses a 
dispersionless Lax pair. We look for a pseudopotential S governed by the equations 

S t = F(S X , u x , u y , u t ), S y = G(S X , u x , u y , ut). (28) 

Calculating the consistency condition St y = S y t and using the equation @ we arrive at five relations 
among F and G, 

Ff-Ga - G^F a = j^F b , 
F^Gb - G^F b = 2j*F b - F a , 

(29) 

F^G C - G^F C = 2^F b + G a , 
G c + ^-F b = 0, G b + 2^F b -F c = 0. 

Job Job 

Here the auxiliary variable £ denotes S x and a = u x , b = u y , c = ut- We also assume f bb ^ 0: for a 
non-trivial Lagrangian, this can always be achieved by an admissible linear change of x, y, t. 

Remark. One has to distinguish between 'true' and 'fake' pseudopotentials, the latter satisfying 
equations of the form 

S t =p(a, S x )c + q(a,S x ), S y = p(a, S x )b + r(a, S x ), (30) 

where the functions p(a,£), q(a,£), r(a, £) solve the system 

Pa + PPi = 0, q a + pq% = 0, r a + pr^ = 0. 

Equations (|3Uj) are automatically consistent by virtue of the relations a y = b x , at = c x , bt = c y 
and, therefore, generate no non-trivial PDEs. Such pseudopotentials can be ruled out, for instance, 
by the requirement F b ^ 0. We assume this hereafter. For a 'true' dispersionless Lax pair (|28|) 
the compatibility conditions must be equivalent to the equation (j2J. One can readily verify that the 
compatibility conditions for the Lax pair (|28jl. (|29|) with F b 7^ imply l|2*j). 

Equations (|29|) imply the expressions for the first derivatives of G in terms of the first derivatives 
of F, 

G c = -&F b , 
c jbb 

G b = F c - 2^F b , 



G * = ~ F ^7Tb ~ F ^ (ft " 2 fc) + 2 fe^ - E t~ 2 fc^< 
G £ = F,(^- 2fyA - 2^ + ^ 



(31) 



along with the following extra relation involving the first derivatives of F only: 

/ cc if if + FlF c (f bb F c - 2f bc F b ) + 2f ac F i F b 2 + 

2F ( F a (f bb F c - f bc F b ) + / bfe F a 2 + f aa F b 2 - 2f ab F b (F^F c + F a ) = 0. 



(32) 



To close this system one proceeds as follows. Differentiating the relation 1)32(1 by a, b, c, £ and calculating 
the compatibility conditions of the equations (|3l|) (that is, Go, = G bc , G ca = G ac , G^ = G^ c , etc, six 
conditions altogether) , one arrives at a linear system of ten equations for the second partial derivatives 
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of F. Solving this system (we skip the resulting expressions for the second order derivatives of F due 
to their complexity) and calculating the resulting compatibility conditions F aa b = F a b a , etc, one arrives 
at exactly the same integrability conditions Q as obtained in Sect 2. from the requirement of the 
existence of hydrodynamic reductions. This proves Theorem 2. 

The necessary details of this computation are explained in the attached file 1-README. The main 
problem was to simplify the expressions for second order derivatives of F: the original expressions found 
using Maple were extremely complicated and would put the calculation of compatibility conditions 
{Faab = Faba-, etc.) beyond our reach: for instance, the number of terms in the expression for F aa was 
equal to 18328, the corresponding number for F ac was 9045, etc. Thus, the hypothetical expression 
for F aa i, would have approximately 10 9 terms. After taking into consideration the identity (|32j) we 
have simplified F aa to an expression with 367 terms only; respectively, F ac had 173 terms. Verifying 
the necessary compatibility conditions and their equivalence to (JUJ) thus became feasible. 

Remark. The functions F and G which determine the Lax pair (|28|) depend on five arbitrary 
constants. Indeed, the values of G, F, F a , -F&, F c can be prescribed arbitrarily at any initial point. On 
the other hand, equations for pseudopotentials remain form-invariant under the transformations 

S — > aS + fix + 7y + 5t + fxu 

which allow one to eliminate this arbitrariness in the generic situation. Hence, for a given non- 
degenerate equation (j2j), the dispersionless Lax pair is essentially unique. 

5 Differential-geometric aspects of the integrability conditions 

Since the integrability conditions Q are preserved by the transformations (el)-(e3), they should be 
expressible in terms of the corresponding invariants. Among the simplest of these invariants is the 
Hessian metric d 2 f. We point out that flat Hessian metrics have been discussed recently in a series of 
publications |IJ [2IJ EH] in the context of equations of associativity of two-dimensional topological field 
theory. Let us recall that, for a Hessian metric d 2 f, Christoffel's symbols and the curvature tensor 
take particularly simple forms, 

Ljj ~f flijt Rijkl P (fimkfpjl fimlfpjk)i 

here f(x l , x n ) is a function of n variables, /y, fyj- denote partial derivatives of /, and the matrix 
p 3 is the inverse of fij. Notice that the Riemann tensor does not contain fourth order derivatives of 
/. The Ricci tensor and the scalar curvature R are defined as R^ = f mp R p i m j, R = f 13 Rji- The 
main object of conformal geometry is the Schouten tensor 

Pi i = ^2 Rij ~ 2(n-l)(n-2) Rflj - 

In three dimensions (n = 3) the condition of conformal flatness is expressible in the form VfcPjj = 
VjPjfc which is equivalent to the vanishing of the so-called Cotton tensor. Our main observation is 
the following 

Proposition. For any density f(a,b,c) satisfying the integrability conditions |0i the Hessian 
metric d 2 f is conformally flat. 

The proof is a straightforward computer calculation. For instance, the density / = abc gener- 
ates the Hessian metric abc (^x + ^l^ + T^) wmcn is manifestly conformally flat. The attached 
program file Cotton. frm performs the general computation in FORM. 

We also present the equivalent tensorial formulation of the integrability conditions ©: 

fijkl = Sym lfijkf Pq fpql ~ ~^fpqifrsjfmkfnl£ P ^ j • (33) 
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Here e^ k is the totally antisymmetric tensor dual to the volume form of the Hessian metric d 2 f (so 
that e 123 = 1/y/H, e 213 = — 1/y/H, etc), and Sym denotes the total symmetrization: 

Sym T ijk i = — T a{i)a{j)a{k)a{l) . 

The formula ()33j) can be rewritten in a completely invariant differential-geometric way as follows. First 
of all, the Hessian metric should be replaced by a metric gij and a flat affine connection V such 
that 

Vfcffij = Vjffjfc, (34) 

indeed, the metric assumes the Hessian form in the flat coordinates of V. Let us define the 3-tensor 
9ijk = ^k9ij (which is automatically totally symmetric) and the antisymmetric tensor e % i k such that 
g i23 _ Yj y/det g^, e 213 = — 1/yMet g~ij, etc. Then (j33j) can be rewritten as 

VfcV; g^ = Sym (gijkg pq g pq i - ^g P qig r sjgm,kgnie prm e qsn ^j . (35) 

Thus, the differential-geometric object underlying the classification of integrable Lagrangians is a pair 
consisting of a metric g^ and a flat affine connection V which satisfy the relations (|34jl. (|35[). 
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